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This  report  compares  calculated  attenuation  rates  for  transverse  magnetic  (TM) 
and  transverse  electric  (TE)  modes  propagating  in  the  earth-ionosphere  waveguide. 
Frequencies  from  10  to  50  kHz  are  considered,  and  attention  is  restricted  to 
situations  where  the  fields  can  be  represented  in  terms  of  the  least- attenuated  TM 
or  TE  mode.  Nine  model  ionospheres  corresponding  to  conditions  caused  by  wide- 
spread high-altitude  fission  debris  are  used  as  inputs  to  the  calculations.  ! 

Certain  of  these  models  are  also  representative  of  polar-cap-absorption  (PCA)  j 
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20.  ABSTRACT  (cont’d) 


events.  The  calculations  utilize  full-wave  methods,  and  account  for  the  verti- 
cal inhomogeneity  of  the  ionosphere,  the  effects  of  heavy  ions,  and  the 
curvature  of  the  earth.  In  addition  to  the  modal  attenuation  rates,  results 
are  given  for  field-strength  height-profiles,  plane-wave  ionospheric  ^flection 
coefficients,  and  the  relative  importance  of  heavy  ions  vie-it-vie  electrons. 
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SUMMARY 


This  report  compares  calculated  attenuation  rates  for  transverse 
magnetic  (TM)  and  transverse  electric  (TE)  modes  propagating  in  the 
earth-ionosphere  waveguide.  Frequencies  fro®  10  to  50  kHz  are  con- 
sidered, and  attention  is  restricted  to  situations  where  the  fields 
can  be  represented  in  terras  of  the  least-attenuated  TM  or  TE  mode. 

Nine  model  ionospheres  corresponding  to  conditions  caused  by  widespread 
high-altitude  fission  debris  are  used  as  inputs  to  the  calculations. 
Certain  of  these  models  are  also  representative  of  polar-cap-absorpt ion 
( PCA ) events.  The  calculations  utilize  full-vave  methods,  and  account 
for  the  vertical  inhomogeneity  of  the  ionosphere,  the  effects  of  heavy 
inns,  and  the  curvature  of  the  earth.  In  addition  to  the  modal  attenu- 
ation rates,  results  are  given  for  field-strength  height-profiles, 
plane-wave  ionospheric  reflection  coefficients,  and  the  relative 
importance  of  heavy  ions  vie-it-via  electrons. 
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PREFACE 

This  report  compares  transverse  electric  (TR)  and  transverse 
nagnetic  (TO)  mode  attenuation  In  the  earth- Ionosphere  waveguide  for 
frequencies  between  10  and  50  kHz.  A wide  variety  of  disturbed  iono- 
spheres characterise ic  of  nuclear  environments  and  polar-cap-absorption 
' PC A ) events  are  considered.  The  results  are  of  interest  in  the  con- 
text of  elevated  antennas  having  a significant  horizontally  oriented 
component . 
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j_.  INTRODUCTION 

I'rvier  most  rood  i t ions , only  transverse  magnet  ic  fTM)  modes  are  -ii 
practical  interest  for  the  propagation  ci  VLF/LF  waves  in  the  earth- 
iorosphere  waveguide.  The  transverse  electric  (Tt)  modes  are  typically 
more  highly  attenuated  and  tore  difficult  to  generate  from  grour.d-fcased 
transmitting  antennas.  ""hus,  fM-morie  attenuation  in  ambient  and  dis- 
turbed environments  has  beer.  thoroughly  studied  and  documented  (e.g., 
.-.nap:  o. . ."-.vaja/'is,  For  elevated,  horizontally 

oriented  transmitting  antennas,  however,  the  TE  mode  could  be  prefer- 
able to  the  TM  mode  under  certain  conditions.  Accordingly,  this  report 
presents  comparisons  between  TF-  and  TX-mode  attenuation  rates  for 
frequencies  between  10  and  50  kHz.  We  consider  a wide  variety  of 
disturbed  ionospheres  characterise ic  of  nuclear  environments  and 
polar-cap-absorptior.  (PCA)  events.  Tn  addition  to  the  calculated 
attenuation  rates,  results  are  given  for  field-strength  height-profiles, 
plane-wave  reflection  coefficients,  and  the  relative  Importance  of  heavy 
ions  ••'a-i-.'fo  electrons. 

Section  II  presents  and  discusses  the  ; or.csphe  r . c models  use.  : 
dec.  LIT,  calculated  TF-  and  TM-moue  attenuation  rates  for  ambient  a: 
disturbed  conditions  is  well  as  f ield-strengt..  height-prof  tie*  for  t:v 
twe  types  of  medes;  Sec-  IV,  TE  and  TM  plane-wave  ret  lei i ion  coefti- 
c ler.ts  for  ambient  and  disturbed  conditions;  and  Sec . 1,  results 
illustrating  the  sensitivity  of  the  propagation  to  ion  ras=,  coll:,.. 
frequency,  and  number  density.  The  equations  and  computer  ro-.es  use-' 
are  summarized  ir  the  appendix. 


I ambient  ant*  disturbed  ionospheric  models 


The  ambient  daytime  electron  and  Ion  density-profiles  used  (see 
Fig.  1)  in  the  calculations  are  taken  from  Knapp  and  Schwartz  (1975). 
The  electron-neutral  collision-frequency  profile  (Fig.  2)  comes  from 
the  same  reference.  The  ion-neutral  collision  frequency,  , is  uncer- 
tain, and  is  believed  to  lie  between  1/10  and  1/40  of  the  electron- 
neutral  collision  frequency,  v , Fxcept  where  otherwise  stated,  the 
calculations  assume  singly  ionized  ions  having  an  atomic  mass  number  of 
25  and  an  ion-to-electron  collision  frequency  ratio  of  1/40. 

The  assumed  nuclear  environments  correspond  to  fission  debris 
uniformly  spread  over  a wide  area  at  altitudes  above,  say,  150  km  or 
so.  The  parameter,  V,  used  to  characterize  the  debris-ionizing  inten- 
sity is  given  by 


W 


A (1+t) 


1.2 


(1) 


where 

Vr  - total  deposited  fission  yield  in  megatons, 

A * area  over  which  debris  is  uniformly  spread 
in  (kilometers) 2 , 

t * time  after  burst  in  seconds. 

Strictly  speaking,  Eq.  (1)  applies  only  to  situations  in  which  all 

bursts  occur  at  t = 0.  For  a large  number  of  bursts  at  different 

1 . 2 

times,  or.e  cart  express  W as  a sum  of  Yp./(l+t)  ’ terms  and  determine  an 
equivalent  value  for  any  time.  However,  the  purpose  of  this  report  is 
to  compare  TK-  and  T>!— mode  propagation  over  a wide  range  of  environments. 
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rather  than  to  draw  detailed  conclusions  for  specific  scenarios.  The 

parameter  V.'  is  a simple  yet  realistic  one  with  which  to  characterize 

—8 

such  a range  of  environments.  For  example,  a value  for  W of  10 
represents  a severe  environment,  values  between  10  ^ and  10  ^ charac- 
terize  moderate  environments,  and  smaller  values  characterize  weak 
environments. 

Figure  3 illustrates  the  normalized,  nominal  ion-pair  production- 

rate  height-profile  due  to  a uniform  layer  of  high-altitude  debris. 

The  peaks  at  about  30  and  65  km  are  due  to  gamma  and  beta  radiation, 

respectively.  The  actual  production  rate  for  a given  environment  can 

14 

be  obtained  by  multiplying  the  curve  in  Fig.  3 by  2.5x10  W;  i.e.,  the 

14  3 

peak  production  is  2.5x10  W ion-pairs/cm  -sec  and  occurs  at  65  km. 

Figure  1 shows  daytime  electron  and  ion  density-profiles  for 
—15  —8 

values  of  !■’  from  2x10  to  2x10  . These  profiles  were  obtained  by 

using  the  production  rates,  as  determined  above,  as  inputs  to  the  quasi- 

equi librium  form  of  the  ionization  balance  equations.  A simple  lumped- 

parameter  model  was  used  for  the  various  reaction-rate  coefficients. 

Although  the  profiles  shown  in  Fig.  1 apply  strictly  to  spread- 

debris  environments,  they  cover  a range  of  ionization  levels  reasonably 

representative  of  other  types  of  nuclear  and  non-nuclear  disturbances . 

-12 

For  example,  the  W = 2x10  profile  is  somewhat  similar  to  that  caused 
by  a strong  PCA  event  (e.g.,  Field , 197G)  . 

Specifically,  F.qs.  (22-25)  to  (22-27)  and  Table  22-5  as  given  I 
All,  Knapp,  and  NUob  are  used  in  the  calculation  of  the  pro- 

files shown  in  Fig.  1. 
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Normalized  Production  Rate 


III.  TE  AND  TM  HOPES 


ATTENUATION  RATES 

The  attenuation  races  of  the  dominant  TE  and  TM  modes  have  been 

calculated  for  each  of  the  model  density-prof lies  shown  in  Fig.  1.  For 

the  more  highly  disturbed  cases,  calculations  are  made  for  frequencies 

between  10  to  50  kHz.  For  the  roughly  2000-  to  10,000-km  path  lengths 

of  interest  here,  the  representation  of  the  field  in  terms  of  only  the 

least-attenuated  mode  becomes  invalid  for  frequencies  above  30  kHz  and 

mildly  disturbed  or  ambient  conditions.  This  lack  of  validity  occurs 

because — unlike  for  moderate  or  severe  environments — in  the  LF  and 

upper  VLF  bands,  the  attenuation  rates  of  higher-order  modes  are  com- 

-14 

parable  with  that  of  the  lowest  mode.  Thus,  for  W £ 2x10  , 30  kHz  is 

the  highest  frequency  for  which  calculations  are  performed. 

Figure  4 shows  the  attenuation  rate  of  the  least-attenuated  TM 
mode  as  a function  of  W for  several  VLF  and  LF  frequencies.  These 
results  correspond  very  closely  to  calculated  TM-mode  attenuation  ratps 
reported  by  CF.  TEMPO  in  Aide  for  the  Study  of  Electromagnetic  Blackout 
(lD7i) . Figure  5 shows  analogous  results  for  the  lea6t-attenuated  TE 
mode.  Comparing  Figs.  4 and  5 reveals  that,  for  all  cases  shown,  the 
TE  mode  is  more  heavily  attenuated  than  the  TM  mode.  In  the  VLF  band 
(<30  kHz),  the  TE-mode  attenuation  becomes  prohibitive  for  environments 
character ized  by  an  ionizing- intensity  parameter,  W,  greater  than  about 
2xl0-11. 

A more  detailed  comparison  between  TM-  and  TE-morie.  attenuation  can 
be  made  from  Figs.  6 through  14,  which  show  the  attenuation  rates  of 
the  lowest-order  TE  and  TM  modes  versus  frequency  for  each  of  the  model 
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profiles  of  Fig.  1.  For  mildly  disturbed  or  undisturbed  conditions 
(Figs.  6 through  10),  the  TE-mode  attenuation  rate  exceeds  that  of  the 
TM  node  by  no  more  than  a few  dS/Mm.  Moreover,  in  the  LF  band  (fre- 
quency L'30  kHs) , the  difference  between  TE-  and  TM-mode  attenuation 
rates  is  only  1 or  2 dB/Mm.  In  these  cases,  the  TE  mode  could  be 
preferred  for  an  elevated  transmitting  antenna  having  a mainly  horizon- 
tal orientation.  For  moderately  and  severely  disturbed  environments 
(Figs.  II  through  14),  the  TE  mcde  is  much  more  heavily  attenuated  than 
the  TM  mode  for  the  frequencies  considered.  In  these  cases,  the  TM 
mode  appears  dominant  except  for  almost  perfectly  horizontal  (electric 
dlpjle)  transmitting  antennas. 

F1ELD-STREMCTK  PROFILES 

To  assess  the  performance  cf  elevated  receivers.  It  is  useful  to 
examine  the  height-profiles  of  the  electric  and  magnetic  field 
strengths.  Accordingly,  these  profiles  are  calculated  for  the  least- 
attenuated  TM  and  TE  modes  for  all  model  ionospheres  and  wave  frequen- 
cies for  which  attenuation  rates  are  shown  above.  In  the  interest  of 
brevity,  only  a few  sample  profiles  will  be  shown  here. 

Figures  IS  and  16  show  the  absolute  values  of  the  calculated 

field-strength  profiles  for  the  TM  and  TE  modes,  respectively. 

-12 

u “ 2x10  and  f ■ 35  kHz  were  chosen  because,  for  these  parameters, 

the  TE-»ode  attenuation  is  only  about  2.5  dB/Mm  larger  than  the  TM-mode 
attenuation  (Fig.  10).  The  r-  and  ^-components  of  the  fields  exhibit 
t'ie  classic  TM-  and  TE-moda  structure,  whereas  the  0-component  exhibits 
a minimum  at  an  intermediate  altitude.  Had  the  real  part  of  the  cal- 
culated fields,  rather  than  the  absolute  values,  been  plotted,  the 
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*1 


4-cor.ponent  height-prc vile  would  exhibit  a more  conventional  structure. 
The  absolute  value  gives  the  aaximum  field  strength  that  could  occur  at 
each  height.  This  taaxiaua,  of  course,  depends  on  the  phase  of  the 
signal  and  will  occur  at  different  phases  for  different  heights. 
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Fig.  4 — Daytime  attenuation  rates  for  the  luwest-order  TM  mode. 


Fig.  ')  -Day tine  attenuation  rates  for  the  lowest-order  TF.  mode 
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Fig.  6 — Attenuation  rate  versus  frequency  for  ambient  day. 
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tenuation  rate  versus  frequency  for  U'  = 2x10 


Atrpmiation  rate  versus  frequency  f nr  V.'  " 2x10 


Fig.  12 — Act 


Fig.  13 — Attenuation  rate  versus  frequency  for  W • 2x10 


Fig.  14 — Attenuation  rate  versus  frequency  for  W - 2x10 
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IV.  REFLECTION  COEFFICIENTS 


Figures  17  through  20  show  the  TM  (E-vector  in  plena  of  incidence) 
and  TE  (E-vector  perpendicular  to  plane  of  incidence)  reflection  coeffi- 
cients at  20  kHz  versus  the  cosine,  C,  of  the  incidence  angle  for 
ambient  conditions  and  disturbed  environments  characterized  by  W - 
2x10  2xl0~^,  and  2x10  respectively.  For  the  ambient  and  mildly 

disturbed  cases  (Figs.  17  and  18),  the  reflection  coefficients  exhibit 
the  classic  behavior;  via,  the  TE  coefficient  decreases  monotonically 
as  C increases,  whereas  the  TM  coefficient  exhibits  a quasi-Brewster 
angle.  However,  each  coefficient  exhibits  a minimum  for  the  more 
strongly  disturbed  environments  (Figs.  19  and  20).  The  reasons  for 
this  anomalous  behavior  are  discussed  below. 

Figures  21  through  24,  which  correspond  to  Figs,  17  through  20, 
show  the  height  dependence  of  the  reflection  coefficients  for  normal 
incidence  (C  - 1)  and  an  oblique  incidence  angle  of  78.46°  (C  - 0.2) 
representative  of  a waveguide  mode  elgenangle.  As  discussed  in  the 
appendix,  j R | is  correctly  interpreted  as  the  ratio  of  dovneoraing  to 
upgoing  waves  only  below  those  altitudes  at  which  reflections  can  occur 
i.e.,  below  the  sensible  ionosphere  for  the  frequency  and  Incidence 
angle  in  question.  Stated  differently,  |r|  is  a true  reflection  coef- 
ficient at  altitudes  where  it  haB  become  independent  of  z.  Thus,  for 
ambient  conditions  (Fig.  71),  significant  reflections  occur  down  to 

_9 

about  55  to  60  km,  whereas  for  W « 2x10  (Fig.  24)  they  occur  at 
altitudes  as  low  as  25  to  30  km. 

To  understand  the  form,  cf  the  reflection  coefficients  shown  in 
rigs.  17  through  20,  it  is  instructive  to  examine  the  height  dependence 
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Fig.  19 — Ret  lection  coefficients  versus  incidence  angle 
angle  at  20  kHz,  W » 2x10“^. 


lertion  coefficients  versus  altitude  for  W ■ 2x10  and  20  kHz. 


Fig.  24 — Reflation  coefficients  versus  altitude  for  W » 2x10  and  20  kHz. 


2 

o*  the  imaginary  part  of  the  squared  refractive  index,  n (see  Eq.  (Al), 

2 

p.  48).  Accordingly,  Fig.  25  shows  /m(n  ) versus  height  for  the  model 
Ionospheres  relevant  to  Figs.  17  through  20.  Field  and  Engel  {1965) 
have  shown  that  most  VLF  reflection  occurs  from  a 5-  or  10-km  wide 
height-region  centered  at  an  altitude  where 

Im  n2  = /2  C2  . (2) 

2 

Moreover,  as  shown  by  Wait  {1970),  if  Im  n has  the  exponential  form 

r 2 z/H  ,,, 

Im  n ^ e , (3) 


v.-here  H is  the  scale-height  of  the  ionospheric  conductivity,  then  the 
TF.-ref lection  coefficient,  Ri,  is  given  by 


-nkHC 


(4) 


vhere  k is  the  free-space  wave  number. 

Ir.  fact,  as  shown  by  Fig.  25,  the  refractive  index  is  locally 

exponential,  having  a scale-height  that  depends  on  altitude.  Consider, 

-9 

for  example,  the  case  W » 2x10  . According  to  Eq.  (2),  the  reflection 

of  a wave  with  C . 0.6  to  0.8  is  affected  by  the  40-  to  50-ka  height- 
range,  where  the  scale-height  is  a maximum.  According  tc  Eq.  (4),  one 
■/ouid  expect  a reflection  minimum  at  these  incidence  angles.  Figure  20 
■■•hitv;:;  that  the  detailed  numerical  calculations  give  such  a minimum, 
conversely,  tor  ambient  conditions,  Eq.  (2)  indicates  that  the  reflec- 
tion altitudes  are  ahove  about  50  km  for  essentially  all  incidence 
angles.  Since  Fig.  25  shows  that  the  ambient  scale-height  is  relatively 
constant  above  50  Itri , no  minimum  would  be  expected  in  the  ambient  day 
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TF-ref lection  coefficient.  The  numerical  results  (Fig.  17)  again  con- 
firm this  heuristic  conclusion.  In  summary,  the  anomalous  behavior  of 
the  reflection  coefficients  under  moderately  or  severely  disturbed 
conditions  occurs  because  waves  with  different  incidence  angles  pene- 
trate to  different  ionospheric  heights,  and  thus  "see"  different 
conductivity  gradients. 

The  results  shown  in  Figs.  17  through  20  must  be  used  with  some 
care  since,  because  of  earth  curvature,  the  incidence  angle  depends  on 
altitude.  In  fact,  since  ionospheric  reflection  does  not  occur  at  a 
well-defined  altitude,  a well-defined  incidence  angle  cannot  be  defined. 
The  best  one  can  do  is  estimate  a height,  h,  that  defines  the  nominal 
center  of  the  reflection  region.  Then  C ae  used  in  Figs.  17  through  20 
refers  to  the  Incidence  angle  at  the  height  h.  (Of  course,  according 
to  Eq.  (2),  h Itself  depends  on  C.)  If  is  the  Incidence  angle  at 
the  ground,  then  the  cosine,  C,  of  the  ionospheric  incidence  angle  is 

C - cos  arc  sinjj—  sin  , (5) 

where  a is  the  earth's  radius;  e.g.,  for  a tangentially  launched  wave 
(«Q  * 90“ ) and  h - 60  km,  C - 0.136. 

It  is  instructive  to  inquire  why  the  TF,  mode  at  20  kHz  is  attenu- 

_9 

atod  by  about  30  dB/Mm  for  W ■ 2x10  , whereas  the  TM  mode  ia  attenuated 
by  only  10  dB/Mm  (Fig.  13).  This  result  may  seem  puzzling  at  first, 
since  lig.  20  shows  that  the  TE-  and  TM-reflection  coefficients  are 
similar  at  oblique  incidence.  A simple  physical  explanation  of  these 
numerical  results  can  be  given  by  interpreting  the  modes  as  plane 
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waves  having  incidence  angles  equal  to  the  real  part  of  the  calculated 


modal  eigenangles. 

To  satisfy  the  boundary  conditions,  the  TE  mode  must  be  more 

-9 

steeply  incident  on  the  ionosphere  than  the  TM  mode.  For  W - 2x10 
and  20  kHz,  for  example,  C was  numerically  calculated  (Eqs.  (A7)  , p.  50, 
and  (A23),  p.  54)  to  be  0.13  and  0.21  for  the  TM  and  TE  modes,  respec- 
tively. One  consequence  of  this  steeper  incidence  is  that  the  TF.- 
reflection  coefficient  is  about  3.5  dB  smaller  than  the  TM  coefficient. 
More  important,  however,  is  the  fact  that  the  steeper  incidence  angle 
causes  the  skip  distance  of  the  TE-mode  plane  wave  to  be  only  about 
350  km,  as  opposed  to  about  680  km  for  that  of  the  TM  mode;*  i.e.,  the 
TE  mode  suffers  tuice  as  many  ionospheric  reflections  per  megameter  as 
+-he  TM  mode.  These  additional  reflections  cause  most  of  the  excess 
TE-mode  attenuation.  Similarly,  much  of  the  anomalous  attenuation  of 
either  mode  in  disturbed  environments  is  caused  by  the  lowered  reflec- 
tion heights  reducing  the  skip  distance.  Thus,  the  attenuation/Mm 
would  increase  even  if  the  reflection  coefficients  were  unaltered. 

For  undisturbed  conditions  or  LF  frequencies,  both  the  TE-  and 
TM-mode  boundary  conditions  are  satisfied  by  very  oblique  launch  angles, 
and  the  TE-mode  skip  distance  is  only  slightly  smaller  than  that  of  the 

TM  mode.  Thus,  for  example,  the  numerical  modal  solu  one  predict  only 

-13 

a moderate  excess  TE-mode  attenuation  for  W - 2x10  and  40  kHz. 


This  skip  distance  is  based  on  a nominal  reflection  height-rani’o 
of  30  to  35  km. 


42 


V.  RELATIVE  IMPORTANCE  OF  IOHS 


For  the  nominal  ion  nan  (29)  end  ion-collision  frequency  (v1/ug 
- 1/40)  used  thus  far  in  this  report,  Eq.  (Al)  (p.  48)  can  be  used  to 
show  that  ion  terms  exceed  the  electron  term  in  the  refractive  index  if 
the  positive  ion  density  exceeds  the  electron  density  by  a factor 
larger  than  700.  Figure  1 (p.  10)  shows  that  this  condition  is  satis- 
fied at  altitudes  lower  than  about  SO  km  for  all  model  profiles  shown. 
Thus,  ions  will  dominate  the  propagation  in  cases  where  the  important 

reflections  occur  below  about  50  km.  Reference  to  Figs.  21  through  24 

-13 

shows  that,  at  oblique  incidence  and  W < 2x10  , the  reflections  occur 

largely  above  SO  km  and  that  iocs  would  be  expected  to  play  a fairly 
minor  role.  Conversely,  for  W Z 2x10  ^ , considerable  oblique  reflec- 
tion occurs  below  SO  km  and  ions  would  be  expected  to  play  a dominant 
role. 

To  quantify  the  above  rather  intuitive  conclusions  regarding  the 
importance  of  ions,  two  types  of  calculations  are  made  with  the  full- 
wave  modal  code  described  in  the  appendix.  First,  the  fraction,  F^, 
of  the  total  joule  heating  in  the  ionosphere  attributable  to  ions  is 
calculated.  Since,  for  the  models  used,  the  attenuation  in  the  wave- 
guide is  due  entirely  to  ionospheric  joule  heating,  represents  the 
fractional  contribution  of  the  ions  to  the  attenuation  rate.  Of  course, 
1 - is  the  fractional  contribution  of  the  electrons.  Second,  the 
modal  attenuation  rates  are  recalculated  using  vi/ve  “ 1/10  rather  than 
1/40.  Actually,  the  results  of  these  calculations  can  be  used  to 
determine  the  sensitivity  of  the  propagation  to  uncertainties  in  ion 
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ic 

mass,  m1,  and  ion  density,  1^,  as  well  as  to  uncertainties  in  v±. 

Since  the  wave  angular  frequency,  u,  is  much  smaller  than  at  the 
altitudes  and  frequencies  of  interest,  the  ionic  contribution  to  the 
refractive  index  (Eq.  (Al))  is  proportional  to  N^/v^m^.  Thus,  a 
f actor-of-f our  increase  in  is  equivalent  to  a factor-of-four  increase 
ic  cij,  a factor-of-four  decrease  in  N^,  factor-of-two  increases  in 
'ind  ra^,  etc. 

Figure  26  shows  Fj  versus  W for  TO  and  TE  modes  and  frequencies  of 

20  and  40  kllz.  As  expected  for  v^/ve  “ 1M0,  ionic  heating  is  the 

-12 

dominant  loss  mechanism  of  W > 10  , and  becomes  minor  only  if  ambient 

-13 

conditions  ate  approached  (W  * 10  ).  For  v^/v^  “ 1/10  (or  equivalent 

changes  in  and  N^) , the  effectB  of  ions  are  considerably  reduced, 
but  still  become  dominant  if  W > 10  11 . Ionic  losses  are  generally 
more  pronounced  for  the  TE  mode  than  for  the  TM  mode,  albeit  by  a 
relatively  small  amount. 

Figure  27  compares  TE-  and  TO-mode  attenuation  rates  at  20  kHz  for 

-12 

;./v  » 1/40  and  v./v  • 1/10.  For  W * 10  , where  electron  heating 

i e i e 

was  seen  from  Fig.  26  to  be  more  important  than  ion  heating,  the 

factor-of-four  increase  in  ion-colll3ion  frequency  produces  only  minor 

changes  in  attenuation  rate.  For  moderate  or  strong  disturbances 
—12 

(w  a 10  ),  the  attenuation  is  considerably  reduced  for  v^/v^  * 1/10. 

Mote  that  the  TF  mode  is  much  more  sensitive  than  the  TO  mode  to  changes 
in  ion  parameters,  with  the  factor-of-tour  increase  in  causing  the 
TE-mode  attenuation  to  change  from  10  dB/Mm  to  5 dB/Mm  for  V » 10  . 

* 

Here,  Ni  is  defined  as  the  combined  number  densities  of  positive 
and  negative  ions. 
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Fig.  26 — Fraction  of  attenuation  caused  by  ion  heating. 


Thus,  the  TE  mode  would  compere  much  more  favorably  with  the  TM  mode  If 
the  iort-coliision  frequency  or  Ion  mess  were  larger  than  assumed  in 
this  report  and  elsewhere  ( Knapp  and  Schxcrtz , 19? S) . Roughly  speaking, 
the  factor-of-four  uncertainty  in  ion-collision  frequency  corresponds 
to  an  order  of  magnitude  uncertainty  in  the  ionizing -intensity  factor,  W. 
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MATHEMATICAL  SUMMARY 

REFRACTIVE  INDEX 

For  computational  purposes,  the  ionosphere  Is  completely  charac- 
terized by  specifying  the  complex  refractive  index  throughout  the 
height-regions  that  govern  long-wave  propagation  in  the  earth- ionosphere 
waveguide.  This  specification  requires  the  nunfcer  densities,  collision 
frequencies,  and  masses  of  each  charged  species  present.  Specifically, 
the  refractive  index,  a,  is  given  by 


2 

n 


1 - 


2 

e 

1837e  m 
o e 


N 

a 

«(w-iva)qa 


<A1> 


where  _ is  the  wave  angular  frequency,  eQ  is  the  electric  permittivity 

of  free  space,  e is  the  electron  charge,  and  m is  the  electron  mass. 

e 

The  number  density,  collision  frequency,  v , and  atomic  mass 
f:  th 

number,  q , of  the  ot  species  can  be  functions  of  altitude.  The 

' £ 

effects  of  the  geomagnetic  field  have  not  been  included  in  Eq.  (Al) 
since,  for  the  disturbed  daytime  environments  considered  in  this 
report,  its  effects  on  mode  structure  are  negligible. 

WAVEGUIDE  MODES 

The  method  used  to  calculate  waveguide  mode  structure  is  based  on 
the  original  formulation  of  Budden  (1961),  and  is  described  in  detail 
for  TM  modes  by  Field  (1970),  Thus,  only  a brief  suanary  will  be  given 


* 

q » 1/1837  for  free  electrons. 
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here.  A spherical-polar  coord inace  eye ten  la  uaad;  the  source  is 
located  at  6 • 0,  and  the  waves  travel  in  the  O-direction.  Because 
azimuthal  (♦)  symmetry  la  assume rl , the  ♦ dependence  is  suppressed 
throughout.  This  assumption  is  valid  because  the  effects  of  the  geo- 
magnetic field  on  the  mode  structure  are  small  for  the  environments 
considered.  Of  course,  if  the  excitation  factor  of  the  nodes  were  to 
be  calculated,  the  $ dependence  would  have  to  be  retained  except  for 

the  case  of  a purely  vertical  electric-dipole  transmitting  antenna. 

1/2 

The  normalized  magnetic  intensity  ff  “ (uQ/eo)  — used,  where  is 
the  magnetic  persdctivity  cf  free  space,  and  B is  the  magnetic  intensity. 
Unless  otherwise  noted,  MKS  units  are  used. 

TM  Modes 

For  TM  modes,  the  electric  and  magnetic  fields  can  be  written 

E(r,«)  - j*rEr<r.e)  + e^r.BjJ  e Ut  (A2) 

and 

jr.O)  -e^(r,e)  eiwt  . (A3) 

For  computational  purposes,  it  is  convenient  to  define  the  wave  admit- 
tance. 


VE« 


(A4> 


and  the  related  quantity. 


(A.-l) 

w*  H Ta^i) 


(A5) 
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By  substituting  Eqs.  (A2)  and  (A3)  into  Maxwell's  equations  and  applying 
the  constitutive  relations,  it  follows  after  some  manipulation  that  W 
is  governed  by  the  following  equation; 


dV, 

dr 


JL  » „2 n u x2 

21  j 


2.2 

! » a-C  ) . 1 

2 2 ,.  2 2 2 

t n 4k  r n 


aw../ 


(A6) 


where  k - u/c,  c is  the  vacuum  speed  of  light,  a is  the  earth's  radius, 
and  n is  given  by  Eq.  (Al).  In  Eq.  (A6) , the  quantity  C can  be  inter- 
preted as  the  complex  cosine  of  the  incidence  angle  of  the  wave  at  the 
ground  (i.e„,  at  r ■ a) . Solution  of  Eq.  (A6)  will  satisfy  the  proper 
boundary  conditions  for  only  discrete  values,  CM,  of  C.  Hence,  C„  is 
called  the  eigenvalue  of  the  TM  mode.  For  a perfectly  conducting 
earth — the  case  considered  here — the  model  equation  for  the  eigenvalue 
is  simply 


W *CN»r"a)  * 1 (A7) 

The  data  given  in  Sec.  II  are  sufficient  to  determine  n“  as  a 

2 

function  of  height.  Once  n is  determined,  Eqs.  (A6)  and  (A7)  form  a 
closed  set  for  W and  C^.,  and  are  solved  by  straightforward  iteration. 
Each  iteration  requires  the  numerical  integration  of  Eq.  (A6) , which  xs 
started  at  a great  height  where  a purely  upgoing  wave  is  assumed  as  an 
initial  solution.  Thus,  since  the  WKB  (elkonai)  solution  may  be  used 
at  preat  heights,  the  assumed  value  of  W is  given  in  terms  of  the 
refractive  index,  n^,  at  the  starting  height,  z , bv 
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(A8) 


W*<*c> 


”o  - V°o  - <1-°2) 

%*  V"o  - (1-c2) 


» 


where  Che  sign  of  the  radical  is  chosen  to  correspond  to  an  upgoing 
wave.  In  practice,  care  is  taken  to  choose  a starting  height  so  large 
that  C is  insensitive  to  its  precise  value. 

Once  CN>  and  hence  Wu(r,CN),  have  been  determined,  it  is  a simple 
matter  to  calculate  all  electromagnetic  parameters  associated  with  a 
given  mode.  First,  it  is  convenient  to  define 


S 


N 


(A9) 


Then,  aside  from  a geometric  spreading  term. 


H, 


can  be  written 


ff^(r-a)  * exp  [ -IkS 


(A10) 


where  D is  the  path  length  along  the  earth's  surface.  It  follows  from 
Fq.  (A10)  that  attenuation  rate,  n,  and  phase  velocity,  V,  are  given  by 


and 


n * 8. 7xl0^k  Im  S 


dB/megameter 


(All) 


V/c  = HRs  S , (A12) 

where  the  .subscript  N has  been  suppressed.  Equations  (AlO)  to  (A12) 
are  valid  beyond  a few  e-folding  distances  from  the  source  and  not  too 
near  the  antipode. 
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The  height-dependences  of  the  fields  are  found  by  a straightforward 
integration  of  Maxwell's  equations,  whence,  suppressing  the  tine  depen- 


dence, 


B (t)  - f-  exp 
<p  r 


11  / 


a JjQ.  dr 

A(r)  dr 


(A13) 


and 


Er(r)  - 


r n fr) 


(A14) 


whereas  is  i.ound  from  Eqs.  (A4)  and  (A5) . The  fields  have  been 

normalized  such  that  B - 1 at  the  ground.  The  aver.  rate  of  power 

♦ 

dissipation  per  unit  volume  is  found  from  Ohm's  lav  be 


V 


P.I  (r)  - -f-  (IeJMeJI)  Im(lV) 


(A15) 


Again,  the  subscript  N has  been  suppressed.  By  using  Eq.  (Al)  in  con- 
junction with  Eq.  (A15),  one  can  determine  the  volume  rate  of  power 
dissipation,  P^Cr),  associated  with  each  of  the  a species.  The  ath 
species  thus  absorbs  a fraction,  F^,  of  the  total  power  lost  by  the 


vth 


mode  to  ionospheric  heating,  where 


F - 

a 


/ V'>  <* 

a 

ow 

f P.t  (r)  dr 


(Ale) 


From  Eqs.  (Al) , (A15) , and  (A16),  it  follows  that 
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TE  Modes 

The  TE  node  fields  are 


E(r,«)  - e E (r,6)  e 
v $ 


iuC 


and 


B(r,Q)  - [erffr(r,e>  + e^r.O)] 


LuC 


Def ine 


and 


- v% 


w*  - 


(At-0 

~ 5 t 

(A;+C‘) 


whence 


dV. 

77 


£ jc^i+wj2  - 4 


2 2 

i - aOkO.  + _ -i 

2 2 ,.2  2 2 

r n 4k  r n 


d-W,  )* 


The  starting  value  for  the  numerical  Integration  of  Eq,  (A21)  is 


w,  (z  ) 

O 


<2  , . „2,i1//2  2 

\no-  U-i  ] j ♦ 

^-u-c2])1'2-,:2 


(A17) 


(A18) 


(A19) 


(A20) 


(A21) 


(A22) 
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The  eigenvalue  of  the  TE  mode  la  found  from 


(Cm,r-a)  - 1 


(A23) 


and  the  formulas  for  attenuation  rate  and  phase  velocity  are  identical 
with  Eqs.  (A9)  through  (A13)  if  is  substituted  for  C^. 

The  fields  are  found  to  be 


(A24) 


B (r ) - - Mr) 

r r 9 


(A23) 


and 


He< r)  - Ax(r)E^(r)  , (A26) 

where  E is  a constant  chosen  to  make  E.  ■ 1 at  the  height  at  which  it 
o $ 

maximizes — typically,  near  the  effective  center  of  the  waveguide  for 
the  lowest-order  TF.  mode.  The  average  rate  of  power  dissipation  is 
given  by 

P,  (r)  - j E~<r)  | Imjl-n2(r)]  . (A27) 


Plane-Wave  Reflection  Coefficients 

The  calculation  of  reflection  coefficient  versus  angle  of  incidence 
suffers  from  interpretive  difficulties  if,  as  has  been  done  above,  the 
full  curvature  of  the  earth  is  included.  In  this  case,  it  is  difficult 
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to  define  an  angle  of  incidence  at  the  ionosphere,  since  the  incidence 
angle  of  a given  wave  changes  with  height  and  the  lower  boundary  of  the 
ionosphere  is  ill-defined.  Thus,  to  permit  a unique  definition  of 
incidence  angle,  the  wave  admittances  that  enter  the  expressions  for 
reflection  coefficients  are  calculated  from  Eqs.  (A5) , (A6) , (A20) , and 
(A21) , but  in  the  limit  at  i + •,  Noting  that  r « a + z,  where  z is 
altitude,  it  follows  that  the  reflection  coefficients  referred  to  z « 0 
are 


C + At(z)  _2ikCz 
R^(2)  - e 


CA29) 


In  these  expressions,  C is  the  cosine  of  the  incidence  angle,  which  is 
to  be  specified  rather  than  calculated  from  the  modal  boundary  condi- 
tions . Strictly  speaking,  R„  and  R;  are  true  reflection  coefficients 
only  at  altitudes  below  those  st  which  reflections  occur;  i.e. , below 
the  sensible  Ionosphere.  At  altitudes  where  the  Inhomogeneous  iono- 
sphere strongly  affects  the  propagation,  the  wave  field  cannot  be 
decomposed  into  downgoing  and  upgoing  plane  waves.  In  practical  terms, 
and  k represent  true  reflection  coefficients  at  altitudes  where 
they  become  independent  of  altitude;  i.e.,  the  mmerical  solutions 
become  independent  of  z (see  Figs.  21  through  24). 
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